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Abstract 

We study the multifield inflationary models where the cosmological perturbation is sourced by light scalar 
fields other than the inflaton. We exploit the operator product expansion and partly the symmetries 
present during the de Sitter epoch to characterize the non-Gaussian four-point correlator in the squeezed 
limit. We point out that the contribution to it from the intrinsic non-Gaussianity of the light fields at 
horizon crossing can be larger than the usually studied contribution arising on superhorizon scales and 
it comes with a different shape. Our findings indicate that particular attention needs to be taken when 
studying the effects of the primordial NG on real observables, such as the clustering of dark matter halos. 
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1 Introduction 

Primordial non-Gaussianity (NG) of the cosmological perturbations has become a crucial aspect of both 
observational predictions of inflationary early universe models and of present and future observational 
probes of the Cosmic Microwave Background (CMB) anisotropies and of the Large Scale Structure (LSS) 
|T|. The main motivation is that detecting, or simply constraining, deviations from a Gaussian distribution 
of the primordial fluctuations generated during an inflationary epoch [2] allows to discriminate among 
different scenarios for the generation of the primordial perturbations. Indeed, a non-vanishing primordial 
NG encodes a wealth of information allowing to break the degeneracy between models that, at the level 
of the power spectra, might result to be indistinguishable. The degeneracy stems from the fact that 
during a period of exponential acceleration with Hubble rate H all scalar fields with a mass smaller than 
H are inevitably quantum-mechanically excited with a final superhorizon flat spectrum. The comoving 
curvature perturbation, which provides the initial conditions for the CMB anisotropies and for the LSS of 



2 



the universe, may be generated not only by the same scalar field driving inflation (the inflaton), but also 
when the isocurvature perturbation, which is associated to the fluctuations of these light scalar fields, 
is converted into curvature perturbation after (or at the end) of inflation [3j|8]. One typical example is 
provided by the so-called inhomogeneous decay rate scenario 6 where the field driving inflation (the 
inflaton) decays perturbatively with a decay rate T. The reheating temperature T r of the hot plasma 
produced by the inflaton decay products is of the order of {M-p\T) l l 2 . If the decay rate depends on some 
light fields which are fluctuating during inflation, then the corresponding large scale spatial variations of 
the decay rate will induce a temperature anisotropy, 5T r /T r ~ ST/F. 

Distinguishing different shapes of the primordial three- (bispectrum) and four-point (trispectrum) 
correlators, i.e. their dependence of the momentum wave- vectors in Fourier space, is of crucial impor- 
tance. Different mechanisms to generate the inflationary perturbations give rise to unique signals with 
specific shapes, which thus probe different aspects of the physics of the primordial universe. For example, 
models in which the curvature perturbation is generated by an initial isocurvature perturbation develop 
(some of the) non-linearities on superhorizon scales. The corresponding NG is of the local type, that 
is the NG part of the primordial curvature perturbation is a local function of the Gaussian part. In 
momentum space, the three-point function arising from such a local NG is dominated by the so-called 
squeezed configuration, where one of the momenta is much smaller than the other two {k\ <C &2 ~ k^). 
The squeezed limit of NG is also particularly interesting from the observationally point of view because 
it leads to pronounced effects on the clustering of dark matter halos and to strongly scale-dependent 
bias [9] . It is impressive that a future detection of a high level of primordial NG in the squeezed configu- 
ration will rule out all standard single-field models of inflation, where the same field drives inflation and 



is responsible for the perturbations, since they all predict very tiny deviations from Gaussianity 10 11 



In Ref. [12] we have taken the first step in trying to characterize the three- and the four-point 
inflationary correlators when the curvature perturbation is generated by scalar fields other than the 
inflaton, the so-called multifield inflation (for the case in which there is only one degree of freedom, see 
Refs. 1 13 1). In particular, we have studied the implications of the symmetries present during a de Sitter 
phase, that is scale invariance and special conformal symmetry. For instance, we have shown that, as 
a consequence of the conformal symmetries, the two-point cross-correlation of the light fields vanish if 
their conformal weights (essentially their masses in units of the Hubble rate) are different. Furthermore, 
we have pointed out that the Operator Product Expansion (OPE) technique is very suitable to analyze 
two interesting limits: the squeezed limit of the three-point correlator and the collapsed limit of the 
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four-point correlator. Despite the fact that the conformal symmetry does not fix the shape of the four- 
point correlators of the light NG fields, we have been able to compute it in the collapsed limit. As we 
mentioned, both the resulting shapes are relevant from the observational point of view. 

In this paper we take a step further and study which informations we can get on the squeezed limit 
of the four-point correlator. This is an interesting question as the four-point function is not fixed by by 
conformal invariance of the de Sitter stage. Once more, we will resort to the OPE technique in order to 
learn what we can say about such a configuration of the four-point correlator. 

One of the goals of these paper is to stress that the contribution to the trispectrum in the squeezed 
limit coming from the NG of the light fields at horizon crossing have a different shape and the amplitude 
can be larger than the trispectrum generated on superhorizon scales (which is there even if the light 
fields are gaussian) . This is somewhat contrary to the common belief spread in the literature whose large 
majority has focused on the NG originated at superhorizon scales. 

The paper is organized as follows. In section 2 we will present a summary of what we know about 
the bispectrum and trispectrum of the light fields during inflation thanks to the symmetry properties of 
de Sitter and the OPE technique. This section contains known material and the expert reader can jump 
directly to sections 3 and 4 where we calculate the four-point correlator in the squeezed limit using various 
arguments. Section 5 contains some quantitative estimates of the three- and four-point correlators from 
the light NG fields; finally section 5 contains our conclusions. 

2 Some general considerations about non-Gaussianities 



We can characterize the cosmological perturbations through the SN formalism 14 , where the comoving 
curvature perturbation ( on a uniform energy density hypersurface at time tf is, on sufficiently large 
scales, equal to the perturbation in the time integral of the local expansion from an initial flat hypersurface 
(t = £*) to the final uniform energy density hypersurface. On sufficiently large scales, the local expansion 
can be approximated quite well by the expansion of the unperturbed Friedmann universe. Hence the 
curvature perturbation at time tf can be expressed in terms of the values of the relevant scalar fields 
ct^ (t* , x) at 



C(i f , x) = Nkj 1 + ^iV/jtr V + ^N IKJ a I a J a K + • • • , (2.1) 
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where Nj, Njj and Njjk are the first, second and third derivative, respectively, of the number of e-folds 

N{t { ,U,x) = [ ' dtH(t,x) (2.2) 



with respect to the field a 1 . From the expansion (2.1) one can read off the n-point correlators. For 



instance, the three- and four-point correlators of the comoving curvature perturbation are given by 



r^MM) = N lNj N K Blf 2%3 , 

B^ n (k u k 2 , k 3 ) = NjNjkNl (P™Pj£ + 2 permutations) (2.3) 



and 

T«-™(k h k 2 , k 3 , k A ) = N lNj N K N L TlfJ^ + N U N K N L N M (pgBffi^ + 11 permutations) , 

T c un (fc!, k 2 , k 3 , k 4 ) = N U N KL N M N N (plLpJMp? N + 11 permutations) 

+ N uk N l N m Nn (p^pJMpKN + 3 permutations) . (2.4) 

Here we have defined kij = {ki + kj) and 



(alal) = (27rf5(k 1 + k 2 )P^ = (2nf5(k 1 + k 2 )6 IJ (^) 

= ( 2 -) 3 ^ + ^ + ^3 + ^Tif^ 4 . (2.5) 

The three-point correlator (and similarly the four-point one) of the comoving curvature perturbation is 
the sum of two pieces 

• One, proportional to the connected three-point correlator of the a 1 fields, is present when the 
fields a 1 are intrinsically NG at horizon crossing. We dub it (in a loose way) the non-universal 
contribution. 

• The second one, which we dub universal (some people come them gravitational), is generated when 
the modes of the fluctuations are super-Hubble and is present even if the a 1 fields are gaussian. 
Of course they depend on the mechanism for the conversion of the isocurvature modes into the 
curvature modes. 
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It is fair to say that most of the attention in the literature has been devoted to the universal contributions. 
For instance, in the large majority of the literature on NG, the nonlinear parameter^ 



fi 



NL 



7NL 



54 

2tnl + t^9nl 
25 



are expressed as 



5 {ChXcX, 



12 pS pS ' 

ki k 2 
1 ^kX^k^kJ 

4 pS p5 pS 



ki k3 ki2 



&4 ' 



P5 ( pS pS +2 permutations 

&4 \ fcl &2 



(squeezed : fci <C &2 ~ ^3), 



(collapsed : &12 — 0), 



(squeezed : <C fci , &2 , k^ 



(2.6) 



5 N^uN- 7 



NL 



6 (iV/iV^) 2 ' 

N I N JI N JK N K 



un 

5nl 



(2.7) 



(iV/iV') 3 ' 
25 N I N J N K N IJK 
54 (NjN 1 ) 3 ' 

with no reference to the contribution to the non-universal terms. This might be due to the fact that the 
non-universal contributions to the connected correlators generated if the light fields are NG depend on 
the specific self-interactions of the light fields and in principle little was known about their magnitude 
and shapes. 

To show that neglecting the non-universal contributions might lead to non correct conclusions about 
NG, let us consider the following simple example where the primordial density perturbations is produced 
just after the end of inflation through the modulated decay scenario when the decay rate of the inflaton is 
a function of some light field a [6], that is T = T(a). If we approximate the inflaton reheating by a sudden 
decay, we may find an analytic estimate of the density perturbation. In the case of modulated reheating, 
the decay occurs on a spatial hypersurface with variable local decay rate and hence local Hubble rate 
H = r(cr). Before the inflaton decay, the oscillating inflaton field has a pressureless equation of state and 
there is no density perturbation. The perturbed expansion reads 



5N d 



Pd 



bn 

3 \p d 



21 



1 The prime denotes correlators without the (27r) 3 5(5^ ki) factors. 
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Immediately after the decay we have radiation and hence the curvature perturbation reads 

<=Wd+ Ht)' (2 - 9) 

Eliminating 5N& and using the local Friedmann equation p ~ H 2 , to determine the local density in terms 
of the local decay rate T = r(<r), we have at the linear order 

c =>(£)• (-0) 

Taylor expanding this expression in powers of the fluctuation a, one obtains 



ffi. = 5 1 



Y'2 



5nl - 54Ar/3 - 3 ^ ^ p/2 + r3 )- (2.11) 

Now, suppose that the function T(a) is of the exponential type, T(a(x,t)) ~ e ao-(^,t) w ^h a constant. 
This is a rather natural possibility if, for instance, the light field is a string modulus setting the amplitude 
of some coupling constant. If so, one immediately concludes that all the universal contributions to the 
NG vanish, = t^l = 5nl = • • • = 0. This holds to any order in perturbation theory as the relation 
between the comoving curvature perturbation £ and the light fluctuation a is linear: £ ~ ln(5r/r) ~ a 



from Eq. (2.10). Alternatively, suppose that the function T(a) gets its dependence on the light field from 
some Yukawa- type interaction with Yukawa coupling Y = Yq(1 + a/2M), with M some high mass scale. 
If so, T(h) ~ r (l + cr/2M) 2 and we obtain negligible NG, = 5/2 and = 25/3, being the current 



bounds / NL ^ 10 2 15 and g NL < 10 6 
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These simple examples indicate that, at least a priori, one may not disregard the non-universal 
contributions unless there is a convincing argument that they are subleading. In fact, it was pointed out 
long ago that the non-universal contributions may be of the same order of magnitude as the universal 
ones [17}{20]. 

From now on we will therefore devote our attention to non-universal contributions to NG in order to 
characterize them. Let us therefore go back to such contributions and briefly summarize what we know 
about them. 

2.1 Non-Gaussianities and the conformal symmetries of de Sitter 

Some progress has been made recently on the knowledge of the NG carried by the light fields during 



the de Sitter stage 112,21,22 . Even though the intrinsically NG contributions to the n-point correlators 



of the light fields are model-dependent, their forms in some specific configurations are dictated by the 
conformal symmetry of the de Sitter stage. Let us recall some of the properties of the conformal symmetry 
in de Sitter. For a more complete description the reader is referred to Ref. [12]. Conformal invariance 
in three-dimensional space M 3 is connected to the symmetry under the group 50(1,4) in the same way 
conformal invariance in a four- dimensional Minkowski spacetime is connected to the 50(2, 4) group. As 
50(1,4) is the isometry group of de Sitter spacetime, a conformal phase during which fluctuations were 
generated could be a de Sitter stage. In such a case, the kinematics is specified by the embedding of 
R 3 as flat sections in de Sitter spacetime. The de Sitter isometry group acts as conformal group on M 3 
when the fluctuations are super-Hubble. It is in this regime that the 50(1,4) isometry of the de Sitter 
background is realized as conformal symmetry of the flat M 3 sections. Correlators are expected to be 
constrained by conformal invariance. All these reasonings apply in the case in which the cosmological 
perturbations are generated by light scalar fields other than the inflaton. Indeed, it is only in such a case 
that correlators inherit all the isometries of de Sitter. First, let us remind how the conformal group acts 
on super-Hubble scales. The set of transformations is given by 

x'i = a i + M i j x j , (2.12) 

x\ = Xxi, (2-13) 

J. - Xi + hi%2 (2 U) 

1 ~ 1 + 2b iXi + Vx* ( ] 

on Euclidean 1R 3 with coordinates x i . These transformations correspond to translations and rotations 

(generated by Pj,Ljj), dilations (generated by D) and special conformal transformations (generated by 

Ki), respectively, acting now on the constant time hypersurfaces of de Sitter spacetime. It should be 

noted that special conformal transformations can be written in terms of inversion 

Xi -> x ■ = ^ (2.15) 

x z 

as (inversion) x (translation) x (inversion). Under conformal transformations, the two-point function of 

fields a 1 and a J of conformal dimensions Aj and Aj respectively, transforms as 

dx' Q x l. Aj/3 , 

{o 1 fry (ft)) -> %± ^ <aV> J (4)> (2-16) 



dxj 



X=X\ 



dXj 



X=X2 



where {dx^/dxjl is the Jacobian of the transformation. For the space inversion (2.15), the two-point 
function, the form of which is forced by scale invariance, transforms as 

{a 1 (x,)a J {x 2 )) -+ 1 Ml 2At ( g f (xi)a J (g 2 )), (2.17) 



2Aj 2A 
x l x 2 



J 



where for x' = x/x 2 we have used that 



dxj 



xfxi 



and the notation x 



i.i 



Xi 



Xj\. Thus, space inversion leaves the two point function invariant if 



Aj = Aj. 

Similarly, the three-point function transforms as 
{a I {x 1 )a J {x 2 )a K {x 3 )) 



dx' { 


A//3 


dx\ 


Aj/3 


dx\ 


dxj 


X=X\ 




X=Xl 





Ax/3 



<aV> J (4)^(4)> 



X=X\ 



and using (2.18), we get that the three-point correlator is invariant if 
A I + Aj - A K , Wl = Aj + A K - Aj 



W K 



wj = A/ + A 



K 



A., 



As a result, two- and three-point function are conformal invariant if they have the form 



(a I (x 1 )a J (x 2 )a K (x 3 )) 



A/ = Aj, 
A/ ^ Aj, 



CIJK 



\x\ - X 2 \ Ai+Aj ~ Ak \x2 - X3\ A J +Ak ~ A '\x3 — xi\ Ai+Ak ~ Aj 



(2.18) 

(2.19) 

(2.20) 

(2.21) 

(2.22) 
(2.23) 



where again a I)J ' K are operators of dimensions Aj jk- In other words, enhancing the symmetry includ- 
ing the special conformal symmetry has two consequences. First, the two-point functions are zero for 
operators with different dimensions and, second, the three-point functions are completely specified by 
special conformal transformations, i.e. by the full conformal symmetry. The four-point function on the 
other hand is not fixed by by conformal invariance. However, as under special conformal symmetry we 
have 



'2 
x 12 



X 



12 



b + xi 



b + X2 



2 ' 



(2.24) 



the four-point correlation can be only a function of the cross-ratios ( t XijX^ rn /xn { Xj rn ). The four-point 
correlator is therefore of the form 



(a {xi)a {x 2 )(T (x 3 )a (x 4 = F , I I x^ 

V Xi 3 X 2 4 ^13^24/ - LJ - J 



-Aj— Aj 



(2.25) 



i<j 



with A = Y^i A/. The four-point function is restricted but not fully specified by conformal invariance to 



be a function of the so-called anharmonic ratios. Therefore, one can conclude that 12 , 22 
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As a consequence of special conformal symmetry the scale invariance during the de Sitter stage, 
the two-point cross-correlation of the light fields vanish if their conformal weights are different. 
Therefore, no assumption is needed on their cross-correlation, it is simply dictated by the conformal 
symmetry [12] . This is the reason why we inserted a Kronecker delta function 5 IJ in the expression 



(2.5). 



The form of the three-point correlator (al ai a*$) is fixed by conformal invariance of the de Sitter 

k\ A?2 &3 

stage and in the squeezed limit it contributes to the total three-point correlator of the comoving 
curvature perturbation (C^C^C^) with the same shape of the universal contributions 



12 



22 



• While the form of the four-point correlator (al ai ai al ) is not fixed by the conformal symme- 

K~L K2 "-3 rC4 

tries of the de Sitter stage, in the collapsed limit it contributes to the total four-point correlator 
^kSk 2 ^k 3 ^kJ w ith the same shape of the universal contributions fl2j: it can be expressed as a 
product of power spectra. 

In the rest of the paper we will characterize the non-universal contributions to the three- and four-point 



correlators expanding the results of Ref. 12 and using the OPE technique. 



2.2 Non-Gaussianities and the Operator product expansion 

The OPE is a very powerful tool to analyze the squeezed limit of the three-point correlator and the 
collapsed and squeezed limit of the four-point correlator. These limits are particularly interesting from 
the observational point of view because they are associated to the local model of NG (for a review see [I] ) 
which leads to pronounced effects of NG on the clustering of dark matter halos and to strongly scale- 
dependent bias (9j. The OPE has been established in perturbative quantum field theories. It is by now 
a standard tool in the analysis of gauge theories such as QCD and Wilson's OPE [23] is the basis of 
virtually all calculations of nonperturbative effects in analytical QCD. It is believed that all quantum 



field theories with well-behaved ultraviolet behavior have an OPE [23 - 25 . 

Let us consider two generic operators cr^(x) and a J (y) at the points x and y on a r = constant 
hypersurface of de Sitter spacetime. Then, we expect that the product of local operators are distances 
small compared to the characteristic length of the system should look like a local operator. As a result, 
we expect that the product of a 1 (x)a J (y) of the two operators cr 7 (x) and a J (y), located at nearby points 
x and y, will have a short-distance expansion of the form [23] 

a\x)a J {y) C n (x - y; g)O n (y), (2.26) 

n 
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where C n {x — y) are c-number functions (in fact distributions), O n local operators and g is the coupling. 
Moreover, for Ht C 1 we expect the OPE to respect the symmetries of the de Sitter spacetime realized 



non-linearly on the r = constant hypersurface. Let us briefly summarize the results obtained in Ref. 12 



Let us note immediately that n-point correlators are reduced to calculation of three-point functions 



by repeated applications of (2.26). Let us consider the case where O is the field itself, i.e. 

a\x)a J {y) ^ £ & K J {x - y; g)a K {y). (2.27) 
K 

The n- and (n + l)-point functions are given by 

g I n 1 ; i In+1 (x 1 ,...,x n+1 -v,g) = (a h (x 1 )---a I ^(x n+1 )y, (2.28) 

9n 



h - I "{x 1 ,...,x n ;n,g) = {a h (x 1 )---a In (x n )Y, (2.29) 



where \x a mass scale. These correlators satisfy the Callan-Symanzik equation 

^^+/3^ + ^7/^/ = 0, (i = n,n + l), (2.30) 

where /3 is the usual /3-function and 7/ the anomalous dimension of a 1 . Using the OPE expansion (2.27) 
one finds immediately that 

9nXi n+1 = E C^gt- 1 ^. (2.31) 
K 

Then, the coefficients of the OPE expansion are also satisfy the Callan-Symanzik equation 

^ + P dg +1I + lJ ~ 1K ) y; /i ' 9) = °' (2-32) 

In particular, for a a conformal field theory for which (3 = 0, dimensional arguments and the fact that 
renormalized operators can be chosen such that they do not depend \x lead to 

CK J (x,y;g) = _ 2wi+2 }j1 2wk , (2-33) 



where u>i t j ; K are the dimensions of the fields a 1 , a J and a K , respectively. Therefore (2.27) can be written 
as 

Here C /l7 should be understood in a non-perturbative sense. In perturbation theory, it has an expansion 

in terms of the coupling(s) g, i.e. for a single field C(g) = Co + c\g + C25 2 -I . Let us now analyze in 

detail the three- and the four-point correlator in various interesting configurations. 
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2.3 The three-point correlator in the squeezed limit 

If we wish to consider the three-point correlator in the squeezed limit, the configuration in real space is 
such that two points, say x\ and x 2 are very close and the third one very far. Let us therefore consider 
the OPE expansion for the two fields a 1 and a J in the (12) channel at the coincident point 

*'(*i)* J (* 2 )= (% + ^*»&) + •••)• (2-35) 
\x 12 x 12 / 



Here w — m 2 /3H 2 -C 1, where m is the mass of the fields, is the conformal weight of the fields involved 
(remember that the weight of the fields a 1 and a 1 must be the same due to the special conformal 
symmetry). The three-point correlator in the squeezed limit can be evaluated by employing the OPE as 

^{x x )a J (x 2 )a K {x 3 )) = ( + ^ A (x 2 ) + • • • ) a K (x 3 )). (2.36) 



Using again the orthogonality of the two-point functions, one finds 12 



K / -> 



C IJ 



A l -At-; \K 



X 



(a A (x 2 )a K (x 3 )) 



12 



C IJK 

x 12 x 23 



(Xi2 C± 0). 



Using the expression 



1 



r( 3=«> 



2«-vr 3 /2r(|) 



d 3 k\k\ 



w—3—ik-x 



(2.37) 



(2.38) 



we obtain for an almost scale invariant spectrum w ~ the Fourier transform of Eq. (2.37) 

(fci < k 2 ~ fc 3 )- 



i a l a l a K\l 



fci k 2 k 3 ' 



QlJKp p 



ki k 2 



1 + ° I % 

ft. 2 



(2.39) 



The non-universal contribution to the three-point correlator in the squeezed limit has therefore the same 
shape of the universal contribution. Its amplitude is model-dependent. 

2.4 The four-point correlator in the collapsed limit 

If we wish to consider the four-point correlator in the collapsed limit, the configuration in real space is 
such that two pairs of points, say x%, x 2 and x 3 , £4 are very far from each other. Let us therefore consider 



the OPE expansion (2.43) as well as the one for the other (34) channel at the coincident point 



a K (x 3 )a^x 4 ): 



C* L (w) , C KL M (w) 



KL , 



„2u> 
"34 



+ 



a M (x 4 ) + - 



'31 



(2.40) 
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The four-point function in the collapsed limit 

{a I (x 1 )a J (x 2 )a K (x 3 )a L (x 4: )) (si2 = and X34 ~ 0) (2.41) 

can be expressed as 



rilj riKL filJ fiKL 

(o*fr)c J &)* K &)o I 'fr)) = -L-ir + i B (a A (x 2 )a B (x,)) + • • • , (2.42) 

x 12 x 34 x 12 x 34 

whose Fourier transforms keeping the connected contribution gives 



^£<«)' ~ ClJ AC KLA P %l P % P %4 + permutations, (k l2 ~ 0). (2.43) 

The non-universal contribution to the four-point correlator in the collapsed limit has therefore the same 
shape of the universal contribution. Its amplitude is model-dependent. 

3 The four-point correlator in the squeezed limit 

Let us now consider the four-point correlator in the squeezed limit which was not analyzed in Ref. |12| . 
We consider three points being close to each other and the fourth very far apart. In other words we 
consider the configuration 

»14 ~ »24 ~ »34 > Vij, (i,3 = 1,2,3). (3.1) 

The method to characterize the four-point correlation in the squeezed limit is based entirely on the OPE. 
As we said, we consider the generic four-point correlator in the squeezed limit in which one of the point 
is much far from the remaining three, say t = x\ 2 — X23 — X13 3> a?i4. We start by dividing space into 
volumes centered around the points x\ , X2 and £3 . For simplicity we will take all of these volumes to 
have the same shape and size R <C £■ We wish to compute the correlator for the smoothed functions 



a 1 R (x i ) = J d'xWRix-x^a^x), (i = 1,2,3), (3.2) 

where Wr(x) is a window function selecting a volume of size R. The key point is that the leading term 
of the four-point correlator in the squeezed limit can be computed considering the modulation of the 
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three-point correlator inside large volume with size at least £ where long wavelength fluctuation modes 
a^(x) live (one can imagine to obtain them by smoothing with a window function (1 — Wr{x)). 
The OPE for the two fields a R and a R in the (12) channel at the coincident point is 

<Wi)<W2) = — + ~a a R ( x v + • • • ' ( 3 - 3 ) 

x 12 X 12 

where we have made explicit that the coefficients of the expansion depend on coupling constants, denoted 
collectively by g and on the long- wavelength mode af(x 2 ) of the field <j n (12). Indeed, this fluctuation is 
seen by the remaining three short wavelength fluctuations as a vacuum expectation value and it should 
be added to the zero mode a N . Multiplying both sides with a R (xs) and since a?i 3 « 0, we get 

4(f!)4(f 2 )af (x 3 ) = ^U(2 3 ) + C " M{W '^ {S2)) + ■■■ (3.4) 

x 12 X 12 

We now expand the three-point correlator C IJ m(w, 9, vf(x 2 )) around a N in powers of erf '{£2) to linear 
order 



4(*iK(* 2 yfe) = g^ogC*,) + ClJMi y ,aN) 



'\2 ^12 



dC IJ M (w,g,a N )aM(x 2 )a«(x 3 )a?(x 2 ) 



da o, 12 
and by applying the OPE in the (23) channel, we get 



<r R {xi)G R {x 2 )(T R (S 3 ) = 2m ^(^3)+ „, 2«, 

x 12 x 12 x 23 

+ is -a <W3) 

12 23 

If we now multiply the above expression by a L (xi) and take the expectation value and Fourier transform, 
we finally obtain 



l„I „J„KL\) 9C IJ m n -\ p /NPK\i, n KL p /_/ _J „M\l /, „ . , . \ /o 7 \ 
^kiVkfkJ ~ ^=i^ C ™ ^ 4 + C ^ P fc 4 ) ' < **> ^ • (3 - 7) 



14 



The non-universal four-point correlator (3.7) can be therefore thought as the modulated three-point 
correlator in the presence of a long wavelength mode at the linear order and at the quadratic order in 
the coupling constants C IJ m 



(al ol a^ah )' ~ l<& (al <ji <£)' A (3.8) 

v fci k 2 k 3 k 4 < \ kx k 2 k 3 ' °7 / v 1 

One can therefore construct first the bispectrum and then replace the (fictitious, if necessary) zero mode 
vacuum expectation value, if any, with a long wavelength mode to be contracted to build up the four- 
point correlator. In particular, the piece linear in the coupling constant C IJ m is similar to the g^L 
contribution, while the quadratic piece in the coupling constant C IJ m is similar to the t-JJl contribution. 
One should remember though that the C IJ m is not only a tree-level quantity, but it is supposed to 
contain informations about all loops. 



The modulation effect in Eq. (3.7) does not come as a surprise. Consider the contribution to the 
four-point correlator of the comoving curvature perturbation coming from the universal terms in Eq. 
(2.4) in the squeezed limit fa <C fa,k2, fa 



T c un (4 fa, fa, k A ) = (2N U N KL N M N N + N IJK N L N M N N ) P£ [P^ ' P? N + 2 permutations) (3.9) 
and the contribution to the three-point correlator of the comoving curvature perturbation coming from 



the universal terms in Eq. (2.3) 



B} n (ki, k 2 , fa) = N^jkNl (Pi K Pi L + 2 permutations) . (3.10) 



Suppose now that there is long wavelength mode, associated to the fluctuation a¥ and to the wavenumber 



A?4, which modulates the bispectrum (3.10). This long wavelength mode is added to the zero mode a M 
from the point of view of the short wavelength modes. One can Taylor expand the derivatives of the 
number of e- folds 



M | M\ 



~ Nj(a M ) + N IM (a M )af, 
* N IJ (o M ) + N IJM {a M )of [ 



(3.11) 



15 



and immediately see that the modulation of the three-point correlator in the background of the long 
wavelength mode gives Eq. (3.9) in the squeezed limit 



T^(k u k 2 , k 3 , h) ~ [agBl* f <7g; fcl, k 2 , k 3 



(3.12) 



Notice that by the same argument one can also reproduce the extra contribution to the non-universal 
trispectrum in the squeezed limit 



rpn—un 



(fci, k 2 , k 3 , k A ) d NijN k N l N m ( Pl« Bd^g^ + 2 permutations 



which comes from the long wavelength expansion (3.11 ) of the term un (fci, k 2 , k 3 



(3.13) 



NjNjNkB 1 /* . 



We would like to stress out at this point that the non-universal contribution to the four-point correlator 



(3.7) is not of the same form of the universal one (3.9). Indeed, in the expression (3.7) the bispectrum is 



not evaluated in the squeezed limit and therefore it is not generically written as the products of power 



spectra. In fact, this is also true for the four-point correlator in the collapsed limit: the expression (2.43) 
is value only in the limit k\ 2 — 0, but is second order in the coupling constants. There might well be 
another piece which does not diverge when k\ 2 ~ 0, but that is first order in the coupling constant. 
Which one dominates clearly depends on the magnitude of such coupling. We will return to this point 
in the next subsection. 



3.1 A consistency check 



To check the validity of the expression (3.8), let us imagine to have only one test field cr(x) with potential 
V{a) = Xa p /p\, with p some positive integer larger than three. Let us also suppose that the light field 
has a vacuum expectation value a which induces an interaction of the form \aa p ~ l / (p — 1)1 (for p = 4 
there will be therefore both a trispectrum and a bispectrum). The n-th correlator point is given by 



K%W- ff « l ( T )> = -''( u 

Using the mode functions in de Sitter 



dr' 



(3.14) 



H 

7W 



(1 — ikr) e 



-ikr 



(3.15) 
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one obtains 
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H 



2n-A ( fc(n) 



n,-2fc 3 



I n (k\, k 2 , • • • , k n ), 



(3.16) 



where ki n) = ki + k 2 + • • • k n , = d n V(o)/do n and 



^cnd 



dr 



Re 



i(l - ifcir) •••(!- ik n T)e lkt 



iki n) r 



(3.17) 



In the squeezed limit k n -C fei, k 2 ,- ■ ■ , k n -\ we have A: n r <C fcjr for (i = 1, . . . , n — 1). Hence, we find that 



In = In-1 + O 



(3.18) 



and the expression (3.16) turns out to be 



R 2n-A j^W 

W2kf~ 



R 2n^ ( k {n-l) 



2k^2kf 



I n [k\, k 2 , • • • , k n ) 

3 

- I n -x(ki, k 2 , ■ ■ ■ , k n -i) 



(or o? • • • a? )' 

\ fcl fc2 fcn-l' 



— -. — — Pj- (or or ■ ■ ■ or ) 

y{n— 1) «n N fcl fc2 fcn-l' 



(3.19) 



In particular, for a potential V(o) = Acr 4 /4!, we find that 

8 Ei<j 2kikj l 



h(h,k 2 ,k 3 ) 
h(h,k 2 ,k 3 ,k4) 



(n) 



(7E + Nfct 



Ei 2Af 



.(3) 



(») 



A; 4 3 



(4) 



3 ' 



(3.20) 



Here 7e is the Euler gamma and Nk t is the number of e-folds from the time the mode kt crosses the 
Hubble radius to the time of end of inflation at T en ^. As a result, the four-point correlator is 

3 

,k 2 ,k 3 ,k 4 ), (3.21) 



H 4 



> fcl ^fc2 *^fc3 ' 

which, in the squeezed limit k 4 <C k\,k 2 , k 3 reduces to 



±U rv-^ n/g 
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("fcl ®ki ® kz ^kA 1 



A 



1 R b% Z-3 Z-3 l3 

±VJ ru-j^ ni/g rt^ 

# 4 ( k\ 



A lfiW k 3 k3 h(h,k 2 ,k 3 ) 



1 



3 ^fei ® k2^ kz I 



k4 

a 



k:i> 



(3.22) 



and reproduces the expression the first piece of the expression (3.8) since in the example at hand we have 



In C IJ ' K = In <7+ constant. It is not difficult to see that the same conclusion holds at second order in 
the coupling constant C (or A) by constructing the bispectrum at this order through the interaction 
Actct 3 and Xa 2 a 2 , see Fig. 1. From the example described above one can also read the information that 



Figure 1: Tree diagrams contributing to the modulation of the three-point function of short 
wavelength modes (continues lines) in the background of a long wavelength mode (dashed lines). 



the four-point correlator in the squeezed limit is not generically written as the products of power spectra, 
as the universal contribution is. Furthermore, at the linear order in the coupling constant A the four-point 
correlator in the collapsed limit does not have the form acquired by the universal piece. This happens 
only at the order 0(X 2 ). Which form is the dominate one depends therefore on the value of the coupling 
A. 

4 The four-point correlator in the squeezed limit: alter- 
native methods 

In this section we wish to offer two alternative methods to get the non-universal four-point correlator in 
the squeezed limit. They are based both on the the symmetries of de Sitter and on the OPE technique. 
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4.1 First method 

The fact that the four-point correlator in the squeezed limit may be obtained from the three-point 



correlator can be also suggested by inspection of the general form (2.25) dictated by the symmetries of 

(4.1) 



de Sitter. Taking the limit of Eq. (3.1), it is easy to see that Eq. (2.25) takes the form 

1 1 pIJKL (*12 X23 \ 

V -''I:'.' -''l.i/ 



Kt=> 



x\X (xi 2 x 23 x 13 ) 2w / 3 
By Fourier transforming the above expression, we get the factorization 



/ a l a l a K a L \/ 

> Jt^-i i*n 



p^ pIJKL 



(k 14 ,k 2 ,k 3 ) 



p^ pIJKL 
fc 4 



(ki,k 2 ,k 3 ), 



(4.2) 



fcl fc 2 &3 &4 

where in the last equality we have used the fact that k\ <C k%, k 2 , k 3 . Then, the function F IJKL depends 
on just three momenta and conformal invariance specifies this function to be necessarily proportional to 



the three point function. In particular, in order to calculate the proportionality factor in (4.2) one can 



make use of the OPEs of the scalar fields. Let us see how this works in a particular model where the 
light fields mix through a quartic interaction of the type 

Mnt ~ XijRLO' I a J a K a L . 



The OPE for the two fields a 1 and a J in the (12) channel at the coincident point is 

oHjB x W&) = + a M (x 2 ) + a\x 2 )a J (x 2 ) + • • • . 

x 12 x 12 

Multiplying both sides with a K (x 3 ) and since x± 3 « 0, we get 

o* = ^IrWa) + ^4^a M (x 2 )a K (x 3 ) 

x 12 X 12 

+ [a I {x 2 )a J {x 2 )y K {x 3 ) + -- - . 
For the operator O ij (x 2 ) = a 1 (x 2 )a J (x 2 ) we have the OPE 



(4.3) 



(4.4) 



O u (x 2 )a K (x 3 ) 



K i 



JjUK ylJKL 



r3w 



+ 



(r L {x 3 ) + 



^23 X 23 

as we assume that the operator IJ a K mixes with a L through the quartic interaction. This is due to in 



2(.c 



(4.5) 



(4.6) 



field theory is called mixing of composite operators. By using Eqs. (4.4), (4.20) in Eq. (4.5) we then get 



I/- \ J(-> \ Ki-> \ M o 

a [xi)o (x 2 )o (x 3 ) = —— 



12 



'12^23 



X 



3uj 
23 



2w 
12 



+ - 



C IJ M C 



MK 



N{W) _jV/-» , X UK N N 



12 



23 



^(*3) + ^-^(*3) 
x 23 x 13 



(4.7) 
(4.8) 
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where for symmetry reasons we have replaced X23 with #13^23 ■ Then, it is easy to find that after 
multiplying the above OPE with a L (x4), the connected part of the four-point correlator turns out to be 

Y)UKL 

{a 1 {x{)a J {x 2 )a K {x z )a L (x^)) = - - 2w + permutations, (4.9) 

X 12 X 23 X 34 

where 

D IJKL ~ X IJKL ~ C U mC MKL I (410) 



By Fourier transforming the expression (4.10) we may get the four-point function in momentum space 



fci k 2 k 3 k 4 < NPM fc 4 \ fej fc 2 £34' 



~ Z^C^P^ <<<^>', (fc 4 « fci,fe,fe). (4.11) 

Therefore, we confirm again that the non- universal contribution to the four-point correlator in the 
squeezed limit does not have the shape of the universal contribution. 

4.2 Second method: Cardy's trick 



For completeness, let us reproduce the expression (3.7) by making use of the conformal symmetries of de 



Sitter. In order to apply the OPE expansion, we need the operators to be at almost coincident points. 
In our case, three points are very close and one point is far form the others. So, to make use of the 
OPE's, we should bring the remote point £4 very close to the points x%,X2 and £3. To do this, we use 



the following trick due to Cardy 26 . We use the fact that under the conformal inversion around an 



arbitrary point xq, the conformal inversion (2.15) becomes 



x' = ,! !V (4.12) 
\x — xo\ 

In other words, if the point £4 is very far from the point £9, then under conformal inversion it becomes 

very close to it. Under the conformal inversion the four-point correlator transforms as 

, j / —* \ 1 / —* \ k,-> •, t,^^ i^ 1 (x\)o- J (x' 9 )a K (x'n)a L (x' d )) ,„„„\ 

(a I (x 1 )a J (x 2 )a K (x 3 )a L (x 4 )) = ^\ 2w )J' } 3 1 ^ , 2w - 4.13 

\X\ — Xq\ zw \X2 — Xq\ ZW \X3 — Xo\ W \X4 — Xq\ zw 

In our case, we may do a conformal inversion about the point xq = £3. Then, the point £3 will remain 
invariant (£ 3 = £3) and the point £4 will come close to £3 as x 34 = I/X34 —> 0. All other distances will 
be clearly larger that x 34 . Thus, we may use the OPE 

zKL nKLMf \ 

a K {x z )a L {^) = + t^-o m {x3) + a K (x 3 )a L (x 3 ) (4.14) 

x 34 ^34 
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where by g we again denote collectively the couplings. We get therefore from the linear term in the above 
OPE the connected part 

C KLM (a I {x' 1 )a J (x' 2 )a M (x 3 )) 



<aV>V 2 )^(*3K(^)) 



X34 \x\ — Xq\ 2w \x 2 — Xq\ 2w \x3 — Xo\ 2w \Xi — Xq\ 2w 



(4.15) 



Now, we go back to the original configuration by doing again a second conformal inversion about the x% 



point. Since the conformal inversion will cancel three terms of the product in Eq. (4.13), we get 



(a I (x 1 )a J (x 2 )a K (x 3 )a L (x i )) 



C 



KLM 



-{a I {x 1 )a J {x 2 )a M {x z )). 



-31 



By Fourier transforming the above expression we finally obtain the second piece of Eq. (3.7) 



(4.16) 



(4.17) 



We should stress again that C should be understood in a non-perturbative sense and have a pertur- 
bative expansion in the coupling constants g, of the form C IJK (g) = ^ n g" ''. 

In fact, it is easy to generalize the above discussion to the general n-point function in the squeezed 
limit. A simple induction of the above leads to 

HA.---°i-siy^ cKLMp kMA~---°i y> (kn«h,k 2 ,...,k n ^). (4.18) 



The first piece of Eq. (|3.7|) emerges from the quadratic term of the OPE in Eq. (4.14). Indeed, this term 

(^1)^2)^3)^3)) 



gives a contribution of the form 

(a I (x" 1 )a J (x' 2 )a K (x 3 )a L (x' 4 )) D 



\£l — Xq\ 2w \x 2 — Xq\ 2w \x3 — Xo| 2u, |x4 — Xq\ 2w 



(4.19) 



For the operator O kl {x 3 ) = a K {x 3 )a L (x^) we have the OPE 
KL {xz)a J {x' 2 ) 



D IJK ylJKL 



\x 2 - x 3 



3w 



x' 2 — X3 1 2w 



(4.20) 



where we assume again that KL a J mixes with a L through the quartic interaction (4.3). Repeating the 
steps above and after an inverse conformal inversion we get 



(^(f 1 ) f T J (x 2 )^(f 3 ) C T L (f4)> D 



ylJKL 
x 12 x 23 x 34 



+ permutations, 



(4.21) 



Therefore the total four-point correlator (3.7) may be written after Fourier transforming as in Eq. (4.11) 



where D IJKL is given by Eq. (4.10). 
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5 On the non-universal contributions to NG 



Having established the form of the non-universal three- and four-point correlator in the squeezed limit, we 
now focus our attention on their magnitude. We already pointed out that the non-universal contributions 
can be the dominant ones. For practical purposes, let us consider the squeezed limit of the three-point 
correlator in a V(a) = Acr 4 /4! model 



- -f i£ P ~A « k * ~ h). (5.i) 

The corresponding contribution to the three-point correlator of the comoving curvature perturbation is 



(Wti' ^-w>ir* Nk < p i p b {kl « k2 ~ h) (5 - 2) 

leading to a non-universal contribution to the nonlinear parameter /nl given by 

fn -un^ 5 _^ N 5 Act N kt _ 5V2X m(a) N kt _ / A \ 1/2 ( m{a)/H \ ( N kt 

~ 1%N'H 2 kt m-KH-pV*- 36vr H pV* ~ \IQ- 2 J \ 1Q- 2 J \ 50 



(5.3) 



where we have defined the quantity = (k 3 /2ir 2 )pS = N' 2 (H/2ir) 2 ~ 2.3- 10~ 9 and the (field dependent) 



k 



Higgs mass m(a) = Act /2. Analogously we find 



5nl - 54-3iY' 2 F2 iVfet " 54 • 3(2tt) 2 P c ~ U \l0- 2 J \ 50 J ' {b ) 

We see that the nonlinearities generated by the non-universal pieces can be substantial, confirming 



previous findings 17 20 . Furthermore, we also deduce that it is quite simple to have models in which 



5nl ^ /nl ( m absolute values), a simple Act 4 model will do it as long as H/a > vA. 

6 Conclusions 

In this paper we have made use of the OPE technique and, partly, of the symmetries of the de Sitter 
epoch, to characterize the NG four-point correlator of the curvature perturbation in multifield inflation. 
In particular we have pointed out that 
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• The contribution to the squeezed limit of the four-point correlator coming from the intrinsic NG of 
the light fields at horizon crossing (which we dubbed non-universal) can be larger than the super- 
horizon contributions (we dubbed them universal) generated even if the light fields are gaussian. 

• The shape of the non-universal contribution to the squeezed limit of the four-point correlator can be 
expressed in terms of the three-point correlator. Nevertheless, in general the shapes of the universal 
and non-universal contributions are different as the three-point correlator is note evaluated in the 
squeezed limit. 

Therefore, particular care needs to be taken when studying the effects of the primordial NG on real 
observables, e.g. the scale dependence of the local halo bias in the presence of NG, as the squeezed limit 
of the four-point correlator is not expressible in terms of products of power spectra. The consequences 
of our findings will be investigated elsewhere. 
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